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Abstract 

A survey is presented of the theoretical status of quadratic response theo- 
ries for the understanding of nonhnear aspects in the interaction of charged 
particles with matter. In the frame of the many-body perturbation theory 
we study the interaction of charged particles with the electron gas, within 
the random-phase approximation (RPA). In particular, nonlinear corrections 
to the stopping power of an electron gas for ions are analyzed, and special 
emphasis is made on the contribution to the stopping power coming from the 
excitation of single and double plasmons. Double plasmon mean free paths of 

swift electrons passing through an electron gas are also discussed. 
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I. INTRODUCTION 



A quantitative description of the interaction of charged particles with matter is of basic 
importance in many different theoretical and applied areasS. When an ion penetrates con- 
densed matter it causes changes in the charge state of the ion, electrons may be stripped 
from the ion or captured from electronic states of the solid, dynamic screening by valence 
electrons originates a wake of electron density fluctuations, and the ion may lose energy 
to the medium through different types of elastic and inelastic collision processes. When a 
swift electron travels in a solid it may also lose energy to the medium. While at relativistic 
velocities radiative losses may become important, for incident charged particles in the non- 
relativistic regime the significant energy losses appear as a consequence of electron-electron 
interactions giving rise to the generation of electron-hole pairs, collective oscillations, and 
inner-shell excitations and ionizations. 

Since the pioneering works of Bohm and Pines! the response of conduction electrons in 
metals to external charged particles has been represented within the electron gas model, by 
replacing the ionic lattice by a homogeneous background which serves to provide neutrality 
to the system. The screening properties of a system of interacting electrons are determined, 
within linear response theory, by the wavevector and frequency dependent longitudinal di- 
electric function e^^^. In the self-consistent field, or random-phase, approximation, the 
dielectric function of an electron gas was first derived by Lindhardi, and, subsequently, a 
number of workers have given alternative expressions for eq,^, incorporating various many- 
body higher order local-field correctionsQ'Q and band effectJi. The effect of dissipative 
processes occurring in a real metal and conversion of plasmons into multiple electron-hole 
pairs may be allowed for in an approximate way by including a damping coefficient in the 
dielectric function!. 

Nevertheless, the validity of linear response theory, which treats the perturbing potential 
to lowest order, is not obvious a priori. Although lowest-order perturbation theory leads 
to energy losses that are proportional to the square of the projectile chargei, Z^e, from 
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measurements on positive and negative pions0 and, also, on protons and antiproton^ll it i 
known that the energy loss exhibits a dependence on the sign of the chargelll~lll. On the 
other hand, experimentally observed nonlinear double plasmon excitations^S cannot be 
described within linear response theoryEBIll, and nonlinearities may also play an important 
role on the electronic wake generated by moving ions in an electron gas00. Finally, lowest 
order perturbation theory breaks down when the projectile is capable of carrying bound 
electrons with itiil. 

The first full nonlinear calculation of the electronic stopping power of an electron gas was 
performed by Echenique et aS, in the low- velocity limit. They used a scattering theory ap- 
proach to the stopping power and the scattering cross sections were calculated for a statically 
screened potential which was determined self-consistently by using density-functional theory. 
These static screening calculations have recently been extended to velocities approaching the 
Fermi velocityii. Alternatively, in the case of incident ions a theoretical effective charge can 
be asso ciatedB, and nonlinearities can be investigated, within quadratic response theory, 
extending, therefore, the range of linear response theory and providing results for arbitrary 
velocity. A quadratic response theory of the energy loss of charged particles in an electron 
gas has recently been carried out@, by following a diagrammatic analysis of many-body 
interactions between a moving charge and the electron gas. 

In this paper we present a survey of the theoretical status of investigations carried out 
within quadratic response theory for the understanding of nonlinear aspects in the interac- 
tion of charged particles with an electron gas. We present general procedures to calculate, 
within many-body perturbation theory, double plasmon excitation probabilities, Zf contri- 
butions to the stopping power of an electron gas for ions and the nonlinear wake potential 
generated by moving ions in an electron gas. We focus on the contribution to the stopping 
power coming from the excitation of single and double plasmons. 

Unless otherwise is stated, atomic units are used throughout {h = nie = = 1) . 
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II. THEORY 



We consider a probe of charge Z\ interacting with a many-particle system. The excitation 
of eigenmodes of the target together with the reaction of the probe to these excitations can 
be described by the self-energy of the probe. For an incoming particle in a state 0o of energy 
p° one writeJni: 



So = I Sy I rfV</.*(r)S(r,r',p°)0o(r'), (2.1) 



where E(r, r',p°) represents the non-local self-energy. 

The real part of Sq gives us the real energy shift due to the interaction with the medium, 
and the imaginary part is well-known to be directly related to the damping rate experienced 
by the particle as a consequence of the interaction with real excitations of the target: 

7 = -2ImSo. (2.2) 

We take the target to be described by an isotropic homogeneous assembly of electrons 
immersed in a uniform background of positive charge and volume Q, and we use, therefore, 
plane waves to describe the incident particle states. Consequently, 

7 = -2ImEp, (2.3) 

where Sp represents the Fourier transform of S(r, r', p"), p = (p,p°), and p is the momentum 
of the probe. 

The self-energy, Sp, can be calculated in the so-called GW approximationlllii: 

^P = ^ZfJ ^,G,.,W„ (2.4) 

where Gk and Wq represent Fourier transforms of the Green function for the probe and the 
time-ordered screened interaction, respectively. In applying this formula we replace Gk by 
the zero order approximation; for electrons {Zi = —1)11: 



where a;k = 77 is a positive infinitesimal, and nk represents the occupation number, 

which at a temperature of T = OK is 

rik = 0{kF - |k|), (2.6) 

kp being the Fermi momentum and 9{x), the Heaviside function. 

The dynamically screened interaction, Wg, can be represented as follows: 

Wg = e-\, (2.7) 

where Vq represents the Fourier transform of the bare Coulomb interaction: 

-. = ^, (2.8) 

and eg is the dielectric function, which is related to the density- density response function, 
Xq, by: 

e-'^l + VqXr (2.9) 

Now, introduction of Eqs. (2.5) and (2.7) into Eq. (2.4), and Eq. (2.4) into Eq. (2.3) 
gives the following result for the damping rate of incident electrons with energy above the 
Fermi level: 



roo 

q 



where Pg represents the probability of transferring four-momentum q = (q, q^) to the electron 
gas: 

47r 

Pq = --^Zf ImWg (5(g° - / + a;p_q)e(u;p_q - Ep). (2.11) 

The delta function in this expression appears as a consequence of energy conservation, and 
the step function, ^(a;p_q — Ep), ensures that no electrons lose enough energy to fall below 
the Fermi level. 
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When the probe is not an electron the occupation number of Eq. (2.6) is zero, i.e., 
we need not take account of the fact that the incident electron cannot make transitions to 
occupied states in the Fermi sea: 

477- 

P, = -—Zl ImW, 6{q' - p° + ^p_q), (2.12) 

and if the probe has mass M » 1, then recoil can be neglected in the argument of the 
delta function to give: 

A'jr 

P, = --^Zllv^WM - q • v), (2.13) 

where v represents the velocity of the incoming particle. 

The inverse mean free path of the probe is easily obtained as follows: 

A-'-iEf^n. (2.14) 

and the stopping power of the target for the probe is obtained as the energy loss per unit 
path length of the projectile, after multiplying the probability Pq by the energy transfer g°: 

In the so-called time-dependent Hartree, or random-phase, approximation the exact lin- 
ear response function to a screened charge is replaced by the response function of the non- 
interacting electron gas: 

XS = -2i / ^.GlGl^q, (2.16) 
thus replacing the linear response function to an external charge, Xq-i by 

xr^=xs+x^.xr^- (2.17) 

Within this approximation the self-energy of Eq. (2.4) can be represented diagrammat- 
ically as in Fig. 1, and cutting the diagrams of this figure through the two-electron lines 
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in all the bubbles would lead to the open diagrammatic representation of scattering ampli- 



tudes shown in Ref. In particular, if M >> 1 the incident particle can be treated as a 
prescribed source of energy and momentum and one findJ^l: 

% = / d^#(g + s - p)W^''^^6{q' - q ■ v), (2.18) 

where s = (s, s°), p = (p,p°), and W^^^ represents the random-phase approximation to the 
screened interaction of Eq. (2.7). 

Then, the probability of transferring four-momentum g to a free-electron gas by moving 
a particle from inside the Fermi sea (|s| < qp) to outside (|p| > qp), thus creating an 
electron- hole pair, is derived from the square of the matrix element S/^i'. 

P, = 2E^sE(l - ^p)l%l'C-., (2-19) 

s p 

where S^^i is the symmetric Kronecker 5 symbol, and introduction of Eq. (2.18) into Eq. 
(2.19) gives exactly the result of Eq. (2.13) found by the self-energy method. 

It is obvious at this point that double plasmon excitations cannot be described within 
the GW-RPA approximation to the self-energy, represented diagrammatically in Fig. 1; 
double excitations can only be described, within the GW approximation, with inclusion in 
the screened interaction of dynamic local-field corrections. On the other hand, the study 
of Zf effects in the stopping power of an electron gas for ions and, also, the study of non- 
linearities in the wake generated by moving ions in an electron gas require going beyond 
the so-called GW approximation. The main ingredient in the investigation of both double 
plasmon excitations and Zf effects is the symmetrized quadratic response function of the 
non-interacting electron gas: 



^k+qi + G^k+q-qi ' (2.20) 



which gives account of the quadratic response of the system to a given charge. The real 
part of this three-point function was first evaluated by Cenni et afl, explicit expressions 
for the imaginary part in terms of a sum over hole and particle states have been presented 
recently0ii, and an extension to imaginary frequencies has also been giveniil. 



A. Double excitation probabilities 



Treating the probe as an external source of energy and momentum, the matrix element 
corresponding to the process of carrying the system from an initial state af^afj^o > to a 
final state o^^aj^l^o > is 



7i/2,ni2 



< *o|a/iO/2 5'a+a+ 


|$o > 


< ^0 





(2.21) 



where $o is the vacuum state, and af are annihilation and creation operators for fermions, 
respectively, and S is the scattering matrix. S is obtained as a time-ordered exponential 
in terms of the perturbing Hamiltonian and field operators ^(x) and '^'^{x) destroying and 
creating, respectively, a particle at the point r at the time t. 

Now, one can apply Wick's theorem, we note that only normal ordered products with 
four uncontracted field operators contribute, and we find, up to second order in the probe 
charge a result that can be represented diagrammatically as in Fig. 2. Within the random- 
phase approximation, the screened interaction, Wg, is obtained from Eqs. (2.7), (2.9) and 
(2.17), and, accordingly, all self-energy and vertex insertions have been neglected. On the 
other hand, exchange processes and, also, ladder contributions have not been introduced 
into Eq. (2.21), since they all lead to scattering probabihties that are of a higher order in 
the screened interaction. 

Finally, the probability for transferring four-momentum g to a free-electron gas by moving 
two particles from inside the Fermi sea (|si| < qp and |s2| < Qf) to outside (|pi| < qf and 
IP2I < Qf) is derived from the square of the matrix element Sf^f^^i^i^: 

Pi = 11(1 - UpJ ^(1 - npj |%/2,iii2l^(5j,p^_,^5^_5^,P2-*2- (2-22) 

qi SI S2 pi P2 

If the probe were not a heavy particle, then recoil should be introduced into the argument 
of the delta function to ensure energy conservation, and, in particular, if the probe were 
an electron an step function should also be introduced to ensure that the probe does not 
lose enough energy to fall below the Fermi level. Then the contribution of Eq. (2.22) 
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to the probability that is proportional to obtained after introduction of the matrix 
element Sf-^f^^i^i^ represented diagrammatically in Fig. 2 into Eq. (2.22), would coincide 
with contributions derived from a GW approximation to the self-energy with inclusion, in 
the screened interaction, of corresponding dynamic local-field corrections. 

In particular, the only Zf contribution to the probability of Eq. (2.22) which might 
represent the real excitation of a double plasmon comes from the square of the scatter- 
ing amplitude represented by the second diagram of Fig. 2. It is given by the following 
expression: 

(2.23) 

Introduction of this probability into Eq. (2.14) gives, after approximating the linear and 
quadratic response functions by their low-g limits, the following high- velocity limit for the Zf 
contribution to the inverse mean free path coming from the excitation of a double plasmonS: 



Numerical study showJll that introduction of the full RPA response functions gives a result 
for which has, in the high- velocity limit, the same dependence on v as the approximation 
of Eq. (2.24), though it is, for = 2.07, larger than this approximation by a factor of 2.16 . 

B. Zf correction to the stopping power for ions 

The stopping power of an electron gas for a probe of charge Zi, mass M » 1 and 

velocity v is obtained after introduction of the probability Pq into Eq. (2.15). Up to third 

order in the projectile charge: 

dE 



I E f ^<l' (n^"^^^ + j ' (2-25) 
where P^^s^*' and p^"^^^'^^ probabilities of transferring four- momentum q to the electron gas 
by creating single and double excitations, respectively, are obtained from Eqs. (2.19) and 
(2.22), respectively. There, 



Sf, = ^Z, j d'qS\q + s - p)S{q' - q • v) 

X {W, + 27rZi [iW,,W,.,,G',+g + -iW,W,,W,.,,Mg,,,] d{q', - qi • v)} , (2.26) 



which is represented diagrammatically in Fig. 3, and 

Sfif2,iii2 = ^^^1 / + Si - Pi) / d^qS\q - gi + S2 - P2)5(g° - q ■ v) 

X { [W,W,^,, + - W,W,,W,.,,M,,,,] + 2mZ,W,,W,.M - Qi • v)} , (2.27) 

represented diagrammatically in Fig. 2. Processes involving higher-order excitations have 
not been included. 

Now, after introduction of Eqs. (2.26) and (2.27) into Eqs. (2.19) and (2.22), respectively, 
we find, up to third order in the ion charge: 

- ^ = (-dEldxr'^^''' + i-dE/dxr'^^"'' + i-dEldxf'''''' , (2.28) 
dx 

where {—dE/dxY^'^^^'^^ represents the contribution to the stopping power coming from 
single excitations: 

i-dE/dxr^'"'^ = --Zl [ r dqy ImW, 6{q' - q • v), (2.29) 

V J (ztt)'^ Jo 

and {—dE/dxY^^^^'^^ and {—dE / dx)'^°'^^^^'^ represent Zf contributions to the stopping power 
coming from single and double excitations, respectively: 

i-dE/dxr^'^^ = --Z! j4\r - q • ^) / ^ /"^ - ^1 • ^) 

V J [ZTr)-^ Jo J [zTr)-^ J-oo 

\mW, Re {W,,W,.,,M,,,,) + Re {w;W,,W,.,,) H,^,,] (2.30) 



X 



and 



ImW„ ImW,_„Im (W,M,„) + Im {w,Wl) ImW,^,, {h,„, + //,„-(,-„))] . (2-31) 



X 



with 
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+ a;s -a;p) - cUg + t^p) 



n^^n n ' (nO n<^] n n + (^l ^ ? " 9l)' (2-32) 

which is related to the imaginary part of the three-point function of Eq. (2.20) bylll'll: 

ImMg^g^ = Hq^q^ + Hq^^q + (2.33) 

The contribution to the Zl stopping power coming from double excitations, which is 
of higher order in the screened interaction than the contribution of Eq. (2.29), has not 
been included in Eq. (2.28). However, contributions of Eqs. (2.30) and (2.31), which are 
proportional to Zf, are all of the same order in the screened interaction, they all need, 
therefore, to be taken into account, and they all can be derived from the knowledge of the 
Zf contribution to the self-energy by going beyond the GW approximation. 

It is interesting to notice that contributions to the stopping power that are propor- 
tional to the product of two imaginary parts of the screened interaction, appearing as a 
consequence of both single and double excitations, can be combined, and we, also, find that 
contributions to the Zf stopping power that are proportional to the product of three imag- 
inary parts of the screened interaction, coming from single and double excitations, cancel 
out. 

Consequently, one finds the following result for the contribution to the stopping power 
that is proportional to Zf: 

X + /2(«,*) + /3(«,«i)l (2.34) 



where 



/i(g, gi) = ImWq ReW^, ReW^.q, ReM,,,,, (2.35) 

/2(g, qi) = ReWq ReWq, ReWq.q, Hq^q,, (2.36) 
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and 



Uq, qi) = -2lmWg ImW,, ReW,.g, H,,,^. (2.37) 

Contributions to the Zf stopping power of Eq. (2.34) coming from /i and /2 of Eqs. 
(2.35) and (2.36) both appear as a consequence of single excitations: the first one comes 
from the cross product between the first and third diagrams of Fig. 3, and gives, therefore, 
the contribution from losses to one-step single excitations generated by the quadratically 
screened ion potential , while the second term comes from the cross product between the 
first and second diagrams of Fig. 3 and gives the contribution from losses to two-step single 
excitations generated by the linearly screened ion potential. The third term comes from 
both cross products, and, also, from losses to double excitations. Contributions coming 
from single plasmons are included in both fi and /s, and contributions coming from the 
excitation of double plasmons are only included in f^. 

Alternatively, the stopping power of an electron gas can be obtained from the knowledge 
of the wake potential induced in the vicinity of the projectile, as the induced retarding force 
that the polarization charge distribution exerts on the projectile itself, and a second-order 
many body perturbation analysis of the wake potentialil at r, defined as the mean value of 
the interaction between a test unit positive charge at that point and the electron gas, leads 
to Eq. (2.34) for the Zf stopping power, as demonstrated in Ref. |25. 

For high velocities of the probe the electron gas can be considered as if it were at rest, 
one can use, therefore, the so-called static electron gas approximation for both linear and 
quadratic response functions, and this results in /2 and /s giving no contribution to the 
integral of Eq. (2.34), i.e., in the high-velocity limit only the contribution to the Zf stopping 
power that is proportional to only-one imaginary part of the linearly screened interaction, 
Wq, is different from zero. Furthermore, it has been shown@ that this contribution to the 
zf effect can be approximated, in the high- velocity limit, by : 

i-dE/dxY^'^ = Zf^Li, (2.38) 
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where Up represents the plasma frequency = Ann, n being the electron density of the 
medium, and Li is the Zf correction to the so-called stopping number: 

Li^l.42^1n— — . 2.39 

At high velocities both the wake potential and the stopping power can also be derived 
within a quantum hydrodynamical model of the electron gas. In this model, we expand the 
nonlinear hydrodynamical equations and find, after quantization, a result for the Hamilto- 
nian of the electron plasma-heavy ion system, in terms of the triple vertex interaction be- 
tween three excitations, that exactly agrees with the result obtained by Ashley and Ritchie0 
by following a different procedure. Then, we find second and third order wake potentials and 
electronic stopping powers, and, also, double plasmon excitation probabilities that coincide 
with plasmon-pole like approximations to the full RPA resultd^. 

III. RESULTS 

Contributions to electron inelastic mean free paths coming from single excitations of the 
electron plasma have been calculated in the high velocity limiti, and, also, in the full RPAii. 
Fig. 4 shows, as a solid line, our full RPA results for the double plasmon inverse mean free 
path of electrons passing through an electron gas of a density equal to that of Aluminum, 
as a function of the velocity, together with the double plasmon inverse mean free path of 
positrons (dashed line) and, also, the high-velocity limit of Eq. (2.24) multiplied by a factor 
of 2.16 (dotted line): 



~ 3.13 X 10-3^. (3.1) 

At high velocities of the projectile the electron gas can be considered to be at rest, the 
effect of the Pauli restriction is, therefore, removed, and the behaviour of the double plasmon 
inverse mean free path, as a function of the velocity, is independent of the particle statistics. 
On the other hand, it is interesting to notice that the high- velocity limit of Eq. (3.1) gives 
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a good account of the full RPA result for both incident electrons and positrons in a wide 
range of projectile velocities. In particular, for Aluminum and an incident electron energy 
of 40keV we find from Eq. (3.1) a ratio for the double relative to the single plasmon of 
1.93 X 10"'^, in agreement with the experiment^. 

Contributions to the stopping power that are proportional to and Zf, as obtained 
from Eqs. (2.29) and (2.34), are plotted in Fig. 5 by solid lines, as a function of the velocity, 
again for = 2.07. It is interesting to notice that both Zl and Zf contributions to the 
stopping power exhibit a linear dependence on the velocity up to velocities approaching 
the stopping maximum; this linear dependence has, also, been observed for the low-velocity 
stopping power when it is calculated to all orders in the probe charge on the basis of density 
functional theoryil. The linear dependence of the Z^ correction to the stopping power is, 
however, a consequence of two competing effects. First, there is the effect of one-step single 
excitations generated by the quadratically screened ion potential, represented by a dashed 
line in the same figure, and, then, the effect of two-step single excitations generated by the 
linearly screened ion potential, represented by a dashed-dotted line. The contribution from 
losses to two-step single excitations, represented by /2 of Eq. (2.36), is very small, at high 
velocities, when the velocity distribution of target electrons can be neglected. In this case 
the static electron gas approximation can be made, the only non-vanishing contribution to 
the zl effect comes, in this approximation, from fi of Eq. (2.35), i. e., from losses to one- 
step single excitations generated by the quadratically screened ion potential, and one finds 
that the result obtained in this approximation is well reproduced by Eq. (2.38), represented 
in Fig. 5 by a dotted line. This approximation gives a good account of the full RPA result, 
even at intermediate velocities where the velocity of target electrons is not negligible, and 
this is, again, a consequence of two competing effects. First, the non-negligible motion of 
the electron gas gives rise to a smaller contribution from losses to one-step single excitations, 
and this is almost compensated by the non-vanishing contribution from losses to two-step 
single excitations. Contributions from losses to double excitations are small in a wide range 
of projectile velocities and they are exactly equal to zero as far as the electron gas can be 
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considered to be at rest. 

In order to analyze the contribution to the nonhnear stopping power coming from losses 
to collective excitations, we first show in Fig. 6 the separate contributions to the linear 
term from plasmon excitation and electron-hole pair excitation by the incident particle. 
For a momentum transfer that is smaller than g^, the critical wave vector where the plas- 
mon dispersion enters the electron-hole pair excitation spectrum, both the plasmon and the 
electron-hole pair excitation contribute to the energy loss, though contributions from losses 
from electron-hole pair excitations are very small. For q > qc, however, only the excitation 
of electron-hole pairs contributes. Total contributions to the Zf stopping power coming 
from q < qc and q > qc are shown in Fig. 7, and contributions coming from q < y/'2ujp and 
q > y/'2u!p are plotted in Fig. 8, y/'2up being the low-density limit of qc- Contributions to 
the stopping power coming from losses to plasmons is, therefore, smaller than contributions 
from losses to electron-hole pairs, especially at high electron-densities, though there is, at 
high velocities, exact equipartition of the energy loss corresponding to momentum transfers 
larger than and smaller than y/2uJ^. This equipartition rule appears straightforwardly in 
the electron gas at rest approximation, and it has been formulated, for an electron gas not 
at rest, by Lindhard and Winther0. This equipartition is, also, found to be exact, in the 
high velocity limit, by using Coulomb scattering of independent electrons with g^m = ^^p/v 
or by assuming that independent electrons are scattered by a velocity dependent Yukawa 
potential with screening length proportional to Up/v. 

As far as the Zf stopping power is concerned, we have split the contributions to /i from 
losses to single plasmons and single electron-hole pairs, and we have found the result shown 
in Fig. 9 by dashed and dotted lines, respectively. On the other hand, all contributions to 
/2, represented in this figure by a dashed-dotted line, come from losses to electron- hole pairs. 
Thus, it is obvious from this figure that contributions to the Zf effect coming from losses to 
plasmons is small, showing that nonlinear corrections to losses from single plasmons are not 
important, and that collective exitations appear to be well described by linearly screened 
ion potentials. The equipartition rule, valid within first order perturbation theory and/or a 
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linear response theory of the electron gas, cannot be extended, therefore, to higher orders 
in the external perturbation. 

Finally, in order to account approximately for the Zf effect coming from both the con- 
duction band and the inner-shells, a local plasma approximation has been used, by assuming 
that a local Fermi energy can be attributed to each element of the solid, and experimental 
differences between the stopping power of silicon for protons and antiprotons have been 
successfully explained in this way§. 



IV. CONCLUSIONS 

In conclusion, we have developed a quadratic response theory for the understanding of 
nonlinear aspects in the interaction of charged particles with matter. In the frame of many- 
body perturbation theory, we have studied the interaction of charged particles with the 
electron gas, within the random-phase-approximation, and in particular, the nonlinear wake 
potential generated by moving ions in matter, the Zf correction to the stopping power, and 
processes involved in multiple excitations of electron-hole pairs and plasmons. 

Double plasmon mean free paths for incident electrons and positrons, and, also, second 
order contributions to the stopping power coming from the excitation of single and double 
plasmons have been evaluated, for the first time, in the full RPA, as a function of the velocity 
of the projectile. 

Our results for the Zf correction to the stopping power show that for velocities smaller 
than the Fermi velocity the stopping power is, up to third order in the ion charge, a linear 
function of the projectile velocity. We have presented, for the high-velocity limit, a formula 
that gives a good account of the full RPA result in a wide range of projectile velocities, 
and our theory gives good agreement with the experiment. We have also separated the 
contributions to the stopping power coming from losses to plasmon generation, and we have 
found that collective excitations are well described by linearly screened ion potentials. 

A nonlinear quantum hydrodynamical model of the electron gas has also been 
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developedEHl. It has been demonstrated that double plasmon excitation probabihties and 
the second order wake potential and stopping power coincide, within this model, with a 
plasmon-pole like approximation to our full RPA scheme, and an extension of this model to 
study the bounded electron gas is now in progres^. 

Full calculations of second order contributions to the wake potential and the induced 
electron density, within the RPA, for different values of the velocity of the projectile and 
the electron density of the medium will be published elsewhere^. 

An analysis of the differences between a self-energy approach to the Zf correction to the 
stopping power and the open diagrammatical approach presented here, and, also, investiga- 
tions of the zf stopping power for incident electrons and positrons are now in progress. 
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FIGURES 



FIG. 1. GW-RPA approximation to the self-energy. 

FIG. 2. Diagrammatic representation, up to second order in the ion charge, of the RPA Sf^f^^i^i^ 
scattering amphtude. Sohd internal lines in the first and second diagrams are zero-order propaga- 
tors, and the triple internal vertex in the second diagram represents the quadratic density response 
function of the non-interacting electron gas. All vertex and self-energy insertions have been ne- 
glected, as well as ladder contributions. 

FIG. 3. Diagrammatic representation of the matrix element of Eq. (2.26), as obtained within 
the RPA. 

FIG. 4. Full RPA double plasmon inverse mean free path of electrons passing through an 
electron gas of a density equal to that of Aluminum (r^ = 2.07), as a function of the velocity (solid 
line). The dashed line represents the double plasmon inverse mean free path of positrons, and the 
dotted line, the high- velocity limit of Eq. (2.24). 

FIG. 5. Full RPA Zi and Zf contributions to the stopping power calculated from Eqs. (2.29) 
and (2.34), respectively, for Zi = 1 and Vg = 2.07, as a function of the velocity of the projectile. 
Dashed and dashed-dotted lines represent Zf contributions from /i and /2 of Eqs. (2.35) and 
(2.36), respectively. The dotted line represents the high- velocity limit of Eq. (2.38). 

FIG. 6. Full RPA zf contribution to the stopping power (solid line), versus velocity. Dashed 
and dotted lines represent contributions from plasmon and single electron-hole pair excitations, 
respectively. 

FIG. 7. Full RPA zf contribution to the stopping power (solid line), versus velocity. Dashed 
and dotted lines represent total contributions for q < Qc and q > qc, respectively, q representing the 
momentum transfer, and qc, the critical momentum for the plasmon being a well-defined excitation. 
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FIG. 8. As in Fig. 7, with qc approximated by its low-density limit: Qc = y^So^. 

FIG. 9. Full RPA Zf contribution to the stopping power (solid line), as a function of the 
velocity of the projectile. The total contribution from /i of Eq. (2.35) (dashed line) has been split 
into contributions coming from losses to single plasmons (dashed-dotted-dotted-dotted line) and 
single electron-hole pairs (dotted line). The dashed-dotted line represents the total contribution 
from /2 of Eq. (2.36), which appears as a consequence of losses to electron-hole pair excitations. 
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